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Abstract 



It has recently been shown theoretically that elastic scattering in the 
Fermi sea produces quantum mechanically entangled states. The mech- 
anism is similar to entanglement by a beam splitter in optics, but a key 
distinction is that the electronic mechanism works even if the source 
is in local thermal equilibrium. An expérimental realization was pro- 
posed using tunneling between two edge channels in a strong magnetic 
field. Here we investigate a low-magnetic field alternative, using mul- 
tiple scattering in a quantum dot. Two pairs of single-channel point 
contacts define a pair of qubits. If the scattering is chaotic, a universal 
statistical description of the entanglement production (quantified by the 
concurrence) is possible. The mean concurrence turns out to be almost 
independent on whether time-reversal symmetry is broken or not. We 
show how the concurrence can be extracted from a Bell inequality us- 
ing low-frequency noise measurements, without requiring the tunneling 
assumption of earlier work. 

To appear in: Fundamental Problems of Mesoscopic Physics: Interac- 
tions and Decoherence, edited by I. V. Lerner et al. (Kluwer, Dordrecht, 
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1. Introduction 

The usual methods for entanglement production rely on interactions 
between the particles and the resulting nonlinearities of their dynam- 
ics. A text book example from optics is parametric down-conversion, 
which produces a polarization-entangled Bell pair at frequency uj out of 
a single photon at frequency 2uj [1]. In condensed matter the schemes 
proposed to entangle électrons make use of the Coulomb interaction or 
the superconducting pairing interaction [2]. 

Photons can be entangled by means of linear optics, using a beam 
splitter, but not if the photon source is in a state of thermal equilibrium 
[3, 4, 5]. Remarkably enough, this optical "no-go theorem" does not 
carry over to électrons: It was discovered recently [6] that single-particle 
elastic scattering can create entanglement in an électron réservoir even 
if it is in local thermal equilibrium. The existence of a Fermi sea permits 
for électrons what is disallowed for photons. The possibility to entangle 
électrons without interactions opens up a range of applications in solid- 
state quantum information processing [7, 8, 9, 10]. 

Any two-channel conductor containing a localized scatterer can be 
used to entangle the outgoing states to the left and right of the scat- 
terer. The particular implementation described in Réf. [6] uses tunnel- 
ing between edge channels in the integer quantum Hall effect. In this 
contribution we analyze an alternative implementation, using scattering 
between point contacts in a quantum dot. We then need to go beyond 
the tunneling assumption of Réf. [6] , since the transmission eigenvalues 
7i,T2 through the quantum dot need not be <C 1. 

The multiple scattering in the quantum dot allows for a statistical 
treatment of the entanglement production, using the methods of quan- 
tum chaos and random-matrix theory [11, 12, 13]. The interplay of quan- 
tum chaos and quantum entanglement has been studied extensively in 
récent years [14, 15, 16, 17, 18, 19, 20], in the context of entanglement 
production by interactions. The interaction- free mechanism studied here 
is a new development. 

The geometry considered is shown in Fig. 1.1. A quantum dot is 
connected at the left and at the right to an électron réservoir. The 
connection is via point contacts connected to single-channel leads. (Spin 
degeneracy of the channels is disregarded for simplicity.) There are two 
leads at the left (Ll, L2) and two leads at the right (RI, R2). A current 
is passed through the quantum dot in response to a voltage différence V 
between the two réservoirs. We consider the entanglement between the 
left and right channels in the energy range eV above the Fermi energy 
Ep. 
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Figure 1.1. Sketch of the quantum dot entangler described in the text. An électron 
leaving the quantum dot at the left or right represents a qubit, because it can be 
in one of two states: it is either in the upper channel (L1,R1) or in the lower chan- 
nel (L2,R2). An example of a maximally entangled Bell pair is the superposition 
(|Ll,Rl) + |L2,R2»/v / 2. 

The degree of entanglement is measured through the violation of a Bell 
inequality [21] for correlators of current fluctuations [22, 23]. Violation 
of the Bell inequality requires mixing of the two outgoing channels at 
each end of the quantum dot (described by 2 x 2 unitary matrices Ul 
and Ur). In order not to modify the degree of entanglement, this inter- 
channel scattering should be local, meaning that it should not lead to 
backscattering into the quantum dot. 1 This might be done by making 
the barrier that séparâtes lead Ll from L2 (and RI from R2) partially 
transparent and tunable by a gâte [23] (cf. Fig. 1.1, shaded rectangles). 

2. Relation between entanglement and 
transmission eigenvalues 

The incoming state, 

i*m)= n a ii°i 2 io), (i.i) 

E F <s<E F +eV 

factorizes into two occupied channels at the left and two empty chan- 
nels at the right, so it is not entangled. Here a} Li (e) is the création 
operator for an incoming excitation at energy e in channel i at the left 
and |0) represents the Fermi sea at zéro température (ail states below 
Ep filled, ail states above Ep empty). There is a corresponding set of 
création operators a Ri at the right. We collect the création operators 
in a vector a) = (a\ 1 , a} L 2 , a R 1: a R 2 ). With this notation we can write 
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the incoming state in the form 

|* m) = at.,.at|0>, „=(Wfv °), (1.2) 

where the product over énergies is implicit. 

Multiple scattering in the quantum dot entangles the outgoing state in 
the two left channels with that in the two right channels. The vector of 
création operators for outgoing states is related to that of the incoming 
states by a unitary 4x4 scattering matrix: b = S- aotf-S = a j! . 
Therefore the outgoing state has the form 

|*out) = tf -S-a-S T -tf\0). (1.3) 



There are two methods to quantity the degree of entanglement of the 
outgoing state: 

A. One can use the entanglement of formation T of the full state 
l'août)- The entanglement of formation of a pure state is defined 
by [24] 

T= -Tr L p L lo g/0L , PL = Tr iï |* out )(* out |, (1.4) 

with Tr£ or Tr^ the trace over the degrees of freedom at the left 
or right. (The logarithm has base 2.) The entanglement of for- 
mation of the outgoing state is given in terms of the transmission 
eigenvalues by [6] 

T = -( e y// i )[r 1 iogT 1 (i-T 2 ) + r 2 io g r 2 (i-T 1 ) 

+ (1 - Ti - T 2 ) log(l - Ti)(l - T 2 )\. (1.5) 

For T\ = T 2 = 1/2 the rate of entanglement production is maximal, 
equal to 2eV/h (bits per second). 

B. Alternatively, one can project |^ ut) onto a state l^out) W1 th a 
single excitation at the left and at the right, and use the concur- 
rence C of this pair of qubits as the measure of entanglement. The 
(normalized) projected state is 

1^/ v = (1 - n L ,in L , 2 )(l - Wfl,iw fl , 2 )|^out) 6) 
' ° ut/ <*out|(l-nL 1 in Li2 )(l-n fl) in fl>2 )|* out ) 1 /2' ^> 

with number operator nx,i = b x { b x i (for X = L, R). The concur- 
rence [25] is a dimensionless number between (no entanglement) 
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and 1 (a fully entangled Bell pair). 2 The transmission formula for 
the concurrence is [6] 

= 2[r 1 (i-r 1 )T 2 (i-r 2 )] 1 / 2 

T l + T 2 - 2TiT 2 ' { ' ' 

Full entanglement is reached when T\ = T 2 , regardless of the value 
of the transmission. 

Notice that in both methods A and B the degree of entanglement dé- 
pends only on the transmission eigenvalues T± , T 2 , and not on the eigen- 
vectors of the transmission matrix. Eqs. (1.6) and (1.7) hold irrespective 
of whether time-reversal symmetry (TRS) is broken by a magnetic field 
or not. In Réf. [6] the expressions were simplified by specializing to 
the tunneling régime Ti,T 2 <C 1. Here we will not make this tunneling 
assumption. 

In what follows we will concentrate on the concurrence C of the pro- 
jected state l^éut)» shrce that is the quantity which is measured by cor- 
relating current fluctuations. The entanglement of formation T of the 
full state |^out) contains also contributions involving a différent number 
of excitations at the left and at the right. Such contributions are not 
measurable with detectors that conserve particle number [26]. 

3. Statistics of the concurrence 

The statistics of C is determined by the statistics of the transmission 
eigenvalues. For chaotic scattering their distribution is given by random- 
matrix theory [11], 

p(t u t 2 ) = cpiTi - r 2 |^(Tir 2 )- 1+ ^ 2 , (1.8) 

with normalization constants c\ = 3/4, c 2 = 6. We obtain the following 
values for the mean and variance of the concurrence in the case (3 = 1 
(preserved TRS) and (3 = 2 (broken TRS): 



^ ~ { 0.3875 if (3 = 2, (L9) 



The effect of broken TRS on the average concurrence is unusually 
small, less than 1%. In contrast, the conductance G = (e 2 /h)Trttf oc 
Ti+T 2 increases by 25% upon breaking TRS. The main effect of breaking 
TRS is to reduce the sample-to-sample fluctuations in the concurrence, 
by about 15% in the root-mean-square value. 
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4. Relation between Bell parameter and 
concurrence 

The Bell parameter £ is defined by [22, 23] 

S = max [E(U L , U R ) + E(U' L , Ur) + E(U L , U' R ) - E(U' L , U' R )] , (1.11) 

where the maximization is over the 2x2 unitary matrices Ul, Ur, U' l , U' r 
that mix the channels at the left and right end of the System. For given 
Ul, Ur the correlator E has the expression 

<(^L,l-^,2)(^fl,l-^fl,2)) ( . m 

((ÔIl,i + ÔI l , 2 )(ÔIr 1+ ÔIr, 2 )}- [ ■ > 

Here SIl,i = Il,ï — is the low-frequency current fluctuation in the 

outgoing channel i at the left 3 and ÔIrj is the same quantity for outgoing 
channel j at the right. The average (• • •) in this équation is over a long 
détection time for a fîxed sample. (We will consider ensemble averages 
later.) 

In the tunneling régime T±,T 2 <C 1 there is a one-to-one relation 
£ = 2\/l + C 2 between the Bell parameter £ and the concurrence C. 
Here we can not make the tunneling assumption. The Bell parameter 
(1.11) can then be larger than expected from the concurrence. The 
relation is [6] 

£ = 2 v / 1 + k 2 C 2 , (1.13) 
k = 1+ - Cft-^) 2 (1.14) 

^Ti(i-ri) + r 2 (i-r 2 ) 1 ; 

The amplification factor k > 1 approaches unity if either T\ « T 2 or 
Ti,T 2 « 1. 

Since £ gives the amplified concurrence kC rather than the bare con- 
currence C, it is of interest to compare the moments of kjC with those of 
C. By averaging with distribution (1.8) we find the mean and variance 
in a chaotic quantum dot: 



{kC) ~ { (L8393 if 13 = % (L15) 
(^)-(nCf = (^38 if = 1, (L16) 



The amplification by k amounts to about a factor of two on average. 
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5. Relation between noise correlator and 
concurrence 

For a différent perspective on the relation between noise and entan- 
glement, we write the correlator (1.12) of current fluctuations in a form 
that exposes the contribution from the concurrence. 

Low-frequency correlators can be calculated with the help of the for- 
mula [27] 

lim {ÔI Ltl {uj)ÔI R , 3 {J)) = -{e*V/h)2itô{uJ+uj')\{rt%\ 2 . (1.17) 

The reflection and transmission matrices r, t are to be evaluated at the 
Fermi energy. We décompose thèse matrices in eigenvectors and eigen- 
values, 

'-^(^ AH t = v «{ y f AH (L18) 

with 2x2 unitary matrices Ul,U r ,Uq. The matrix r contains the re- 
flection amplitudes from left to left and the matrix t contains the trans- 
mission amplitudes from left to right. 4 
Substitution into Eq. (1.12) gives 

E = (l-2\U L ,n\ 2 )(l-2\U R ,n\ 2 ) + 4 K CReU L<n U* Rjn Ul 12 U RA2 . (1.19) 

We see that the entire dependence of the correlator E on the transmission 
eigenvalues is through the product kC of concurrence and amplification 
factor. This is the same quantity that enters in the Bell parameter 
(1.13). The correlator E is less useful for the détection of entanglement 
than the Bell parameter S, because it dépends also on the matrices of 
eigenvectors Ul,U r — which the Bell parameter does not. 

In a chaotic quantum dot the two matrices Ul and U R are indepen- 
dently distributed in the circular unitary ensemble (a socalled "isotropic" 
distribution [11]). Averages over thèse matrices can be done conveniently 
in the parametrization 

jj _ e ia l e ^ cos 7 e ^ sm 7 
l _g-# s i n ^ e _î< ^cos7 

7 G (0,7r/2), a,<j),ij£ (0,2vr). 

The isotropic distribution implies that ail four angles 7, a, (f>, ip are inde- 
pendent. The distribution of a, <f>, ip is uniform while the distribution of 
7 is P(j) oc sin27. 

With this parametrization Eq. (1.19) takes the form 

E = cos2j L cos2j R + KCsm2j L sin2-f R cos((pL-tpL-<pR + ipR)- (1-22) 



(1.20) 
(1.21) 
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Figure 1.2. The quantum dot of Fig. 1.1 has been replacée! by a disordered wire 
(dotted rectangle). Although the distribution of transmission eigenvalues is différent, 
the relation (1.23) between noise correlator and concurrence still applies. This relation 
only relies on the isotropy of the eigenvector distribution. 

Upon averaging over the angles we find 

(E) = 0, (E 2 ) = \ + \(k 2 C 2 ). (1.23) 

The significance of this équation is that it applies generally to 2 x 2 
transmission matrices with an isotropic distribution of eigenvectors, even 
if the distribution of eigenvalues differs from Eq. (1.8). For example, it 
applies to the disordered conductor shown in Fig. 1.2. 

6. Bell inequality without tunneling assumption 

The Bell parameter £ is no longer directly related to the concurrence 
C if the transmission probabilities are not small compared to unity [6]: 
The relation (1.13) between £ and C contains a spurious amplification 
factor k > 1, which approaches unity in the tunneling régime. The same 
amplification factor appears in the correlator (1.19). In this section we 
show how one can avoid the amplification factor, by calculating the vio- 
lation of the Bell inequality without making the tunneling assumption. 
The same problem was studied, from a différent perspective, in Réf. [9]. 

As described in Réf. [22], the Bell inequality is formulated in terms 
of the correlator Kij of the number of outgoing électrons detected in a 
time r in channel i at the left and channel j at the right: 

= r- 2 r dt r dt' (i L , t (t)i R>J (t')) 

Jo Jo 

^du ^ l ' C^co). (1.24) 
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Here Cij{uj) is the frequency dépendent correlator of current fluctuations, 

/oo 
dté^(ÔI Lti {t)5I R ^)). (1.25) 
-oo 

In the tunneling limit it is possible to neglect the product of averages 
(lL,i)(lR,j) and retain only the second terni in Eq. (1.24), proportional 
to the current correlator dj . Both terms are needed if one is not in the 
tunneling limit. 

We assume that V is small enough that the energy dependence of the 
scattering matrix may be neglected in the range (Ep,Ep + eV). (That 
requires eV small compared to the mean level spacing of the quantum 
dot.) Then the frequency dependence of Cy(cj) is given simply by 5 

For short détection times r <C h/eV one may take the limit 



In view of Eq. (1.17), the zéro- frequency limit of the current correlator 
(1.25) is given by 



C t3 (0)=-(e 3 V/h)\(rt%\ 



(1.28) 



The mean outgoing currents are given by {Il,ï) = (e 2 V '/h)(rr^)a and 
(Irj) = (e 2 V/h)(ttî)jj. Substitution into Eq. (1.24) gives the short- 
detection-time limit 



limtftf = {I Lyl ){I R ,) + {eV/h)C^) 



(e 2 V/hf \(rr^Utt%-\(rt%\ 



We now define the correlator E in terms of the short-time Kij, 

K u + K 22 - K 12 - K 21 



E = 



K n + K 22 + K 12 + K 21 ' 



(1.29) 



(1.30) 



Notice that this définition of É corresponds to définition (1.12) of E if 
is replaced by CV,(0). Substitution of Eq. (1.29) leads to 



^ = -(l-2|C/ Lill | 2 )(l-2|^, 11 | 2 )-4CReC/ Lill f/^ 11 C/£ il2 [/ Ril2 , (1.31) 
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where we have used the parametrization (1.18). Apart from an overall 
minus sign, Eq. (1.31) is the same as Eq. (1.19) — but without the factor 
k multiplying the concurrence. 

The maximal violation £ of the Bell inequality is defined in the same 
way as in Eq. (1.11), with E replaced by E. The resuit 

£ = 2y/l + C 2 (1.32) 

is the same as Eq. (1.13) — but now without the factor k. 

Since short-time détection experiments are very difficult in the solid 
state, the usefulness of Eq. (1.32) is that it allows one to détermine the 
concurrence using only low-frequency measurements. It gêner alizés the 
resuit of Réf. [23] to Systems that are not in the tunneling régime and 
solves a problem posed in Réf. [6] (footnote 24). 

7. Conclusion 

We have investigated theoretically the production and détection of 
entanglement by single-electron chaotic scattering. Much is similar to 
the tunneling régime studied earlier [6], but there are some interesting 
new aspects: 

■ The degree of entanglement, quantified by the concurrence C, is 
sample spécifie. The sample-to-sample fluctuations become smaller 
if time-reversal symmetry is broken, while the average concurrence 
is almost unchanged. 

■ The low-frequency current correlator Cy and the Bell parameter 
£ constructed from it give the concurrence times an amplification 
factor k. In the tunneling régime n 1. One also has n = 1 
if the two transmission eigenvalues Ti,T% are equal. The factor 
k can become arbitrarily large if T\ — ► 1 and T2 — > (or vice 
versa). On average, the amplification factor in an ensemble of 
chaotic quantum dots is about a factor of two. 

■ The bare concurrence, without the amplification factor, is obtained 
by adding to the low-frequency current correlator the product of 
average currents times h/eV . 

■ The concurrence gives the maximal violation of the Bell inequality 
for détection times r short compared to the cohérence time h/eV. 
In Réf. [22] the opposite limit r S> h/eV was taken, which is 
appropriate in the tunneling régime, but does not allow to violate 
the Bell inequality outside of that régime. A similar conclusion 
was reached in Réf. [9]. 
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From an expérimental point of view, the missing building block in Fig. 
1.1 is the local mixer at the left and right end of the quantum dot. Thèse 
mixers are needed to isolate the contribution to the noise correlator that 
is due to the concurrence. In optics, a simple rotation of the polarizer 
suffices. The electronic analogue is a major challenge. 
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Notes 

1. The mixers have no efïect on the incoming state, because both incoming channels are 
either filled or empty at any given energy. 

2. The concurrence C of the qubit pair is related to the entanglement of formation T 1 of 
the projected state |^o ut ) by T' = —xlogx — (1 — x)log(l — x) with x = i + ^ Vl — C 1 . 

3. The total current in channcl i at the left (incoming minus outgoing) is e 2 V/h — lL,i- 

4. In the présence of TRS one has Uo = Uj^ , but this constraint is irrelevant because 
anyway Uo drops out of Eq. (1.17). 

5. The cross-correlator C'ij(u) vanishes for \hui\ > eV because wc are correlating only the 
outgoing currents; the correlator of incoming plus outgoing currents contains also a voltage- 
indcpcndcnt tcrm oc \hu>\, cf. Réf. [28]. 
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